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ABSTRACT 


Vertical stern accelerations for five geometri-~ 
cally similar destroyer hull forms of different lengths 
in the 250-550 foot range are predicted on the basis of 
model tests for a ship speed of 17 knots at various 
headings in representative shortecrested irregular seas 
of varying severity. 

Using tne concept of statistical inference, the 
model test predictions are treated as a statistical 
sample and used to predict what would occur in the total 
population of possible stern accelerations. 

Then, with strong emphasis on the systems 
approach, it is shown how these predictions can be used 
as input to a probability analysis thet leads to a pre- 
Giction of long range system behavior/effectiveness, 
the system under consideration being the Destroyer/ 


Variable Depth Sonar System. 





DD 


RAO 


7) 


Effective7 
i 


Amplitude 


Height 


el 


oe ot 


byl 


NOMENCLATURE 


Destroyer type shin. 
Variable Depth Sonar. 
Response Amplitude Operator. 


Actual wave length of regular waves in 
feet measured perpendicular to lines of 
eresgus . 

Wave length apparent to a ship proceeding 
obliquely through regular waves. 


Ship length in feet measured along 
designer's waterline 


For a regular or irregular function, 
peak to peak value civided by two. 
Sometimes referred to as single 
amplitude. 

For a ré@ular or irregular function, 
peak to peak value. Sometires referred 
to as double amplitude. 

Amplitude of ship pitch angle in radians; 
vertex of angle is defined as being at 
ship's center of gravity. 

Instantaneous value of pitch. 

Amplitude of heave (vertical transla- 
tion) of shin in feet at center of 
gravity. 

Instantaneous value of heave. 

Amidships 


Vertical translation at a point x feet 
from B. 


Wave amplitude in feet. 


Li 





can a 


me 2,712 


Instantaneous value of distance in feet 
between wave surface and calm water zero 
datume 

Maximum wave slope. see Page 9. 


Phase angle between pitch and wave, 
See Pages 8 and 9. 


Phase angle between heave and wave. 
see Pages 8 and 9. 


Phase angle between heave and pitch. 
See Page 9. 


See Figure 3 on page 11 for angular 
nomenclature, 


Ship speed in feet per second, 


Circular frequency of wave in radians per 
second. 


Encounter frequency in radians per second, 
Time in seconds, 


eee cee on due to gravity; 32.2 ft. per 
(seconde ) 


Distance along ship's length from & ~~ 
some point, (+) forward, le ) eft. 
Appendix A, ordinate of irregular oe 
acceleration record, 


Amplitude of vertical acceleration at a 
point x along ship's length. In this 
thesis = -181 feet for the 370 foot 
ship length. (x changes with geometric 
similarity for other ship lengths) 
Logarithm to the natural base e, 
Increment of. 


Piers on-Moskowitz sea spectrum ordinate, 
"Corr" implies corrected. See Page 16 


Lag number; sea spectrum abscissa fre- 
quency coordinate. See Pages 16 and 17, 
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Sample 
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N(x,8”) 


Population 


Yo f 
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Significant wave height in feet; defined 
as the average of the 1/3 highest wave 
heights extant in any given seaway. 


Response Spectrum Component. See Figure 9. 
Integrated Response Spectrum. See Figure 9. 
Ship response parameter. See Appendix A. 


A relatively small portion of the set of 
values in the universe of a random variable. 


Mean of sample probability distribution of 
stern vertical acceleration. 


Variance of sample probability distribution 
of stern vertical acceleration. 


Deviation of sample probability distribu- 
tion of stern vertical VE acceleration. 


A normal (Gaussian) probability distribu- 


tion function whose mean is X and whose 
variance is S°“. 


The set of all possible values in the 
universe of a random variable, 


Mean of population probability distribution 
of stern vertical acceleration. 


Variance of population probability distri- 
bution of stern vertical UE acceleration; 
when used in a general sense, the variance 
of any normal probability distribution 
function. 


Deviation of population probability distri- 
bution of stern vertical VE acceleration. 


True population parameters. 

Probability of [ ]. 

Estimate of “o, So. See Figure 17 on Page 56. 
Location of X and s in S/W parameter space. 


Size of a statistical sample. 
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R z 
RMS * 
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Suverscripts 


nn 
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ea 
Subscripts 

Q - 
max - 
- om 
y os 
mT ~ 
>: @ a) 
e - 
Ba ~ 
ig é 
PTP ~ 


Number of cycles of a certain ship motion. 


Parameter of the Rayleigh probability dis- 
tribution function. 


Root mean square of some variable. See 
Appendix A. 


See Aopendix A. 


Denotes amplitude of a sinusoidal 

function as distinguished from instan- 
taneous value. In the case of an irregular 
function, denotes peak to peak divided by 
two. This superscript does not imply the 
above when used in X; in that case it 
indicates an arithmetic mean. 


Average of. 


Measured at shin's center of gravity. 
Maximum value. 

heave relative to wave. 

pitch relative to wave. 

pitch relative to heave. 


at location x feet along ships length 
from 8. 


Encounter. 

the {th value of a variable. 
upper. 

lower. 


peak-to-peak value (double amplitude or 
height ) 





the highest 1/3 of a group of values 
of a variable. 


original value, before being changed. 


value after having been changed from 
original value. 
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INTRODUCTION 


The basic method used in this thesis to predict 
the various modes of response of a particular ship to a 
given seaway was introduced over ten years ago (1)* and 
has been extensively employed in the intervening time. 
The challenge of the present day is how to most effec- 
tively utilize such response data in order that the 
systems designer might be provided with reliable long- 
range estimates of system behavior and as well, that the 
system user might have reliable estimates of system 
effectiveness, 

The approach developed in this thesis is a 
broad and general one that is applied to one particular 
aspect of behavior of one system. It is suggested that 
the reader keep in mind the generality of the method, 
while considering the specific case which is developed 
herein for illustration as well as for its current 
naval interest, 

The system considered in this thesis is the 
DD/VDS system. The particular aspect of system behavior 
considered is the vertical acceleration of the VDS tow- 


point at the destroyer's stern. 





“Numbers in parentheses refer to references listed on 
Page 91. 





The development undertaken herein can be divided 


into three principal areas, namely: 


(1) prediction of vertical stern acceleration 
(in terms of the VE parameter which is 
discussed in detail in Appendix A), for 
various ship lengths, headings and sea 
States; 

(2) utilization of the concept of statistical 
inference to predict the population dis- 
tribution of vertical stern acceleration™ 
(in terms of the VE parameter), treating 


the predicted values in (1) above as a 





statistical sample; 

(3) a demonstration of how these predicted stern 
accelerations can be used in a probability- 
based system analysis leading to long- 
range prediction of DD/VDS system behavior 


over the lifetime of such a system. 





“Phe term "population™ is used here to signify all the 
possible values of vertical stern acceleration (in terms 
of the VE parameter) that a destroyer of the particular 
type considered might experience in its lifetime. In 
theoretical probability and statistics, “population” and 
"universe" are used interchangeably to describe the set 
of all the possible values of a given random variable, 








PREDICTION OF STERN ACCELERATION RESPONSE 


MR eS ee 


General Procedure for Prediction of Sina Re sponse 
In 1953, St. Denis and Pierson (1) introduced a 


method utilizing the principle of superposition whereby 
a ship's response in irregular waves may be preaicted 
provided the response in regular waves is known, as it 
is, for examole, from model tests in regular waves. 

The method basically consists of (a) representing 
the sea as an energy spectrum, (b) transforming model 
response data taken in regular waves into a response 
amplitude operator (RAO) (referred to by some as a 
transfer function) and (c) multiplying the RAO by the 
Bea spectrum to arrive at a@ response spectrum having 
certain probability properties regarding the response 
being consicered. 

One section each will be devoted to the descrip- 
tion of develoving RAO's from model data, to the selection 
of sea spectra, ana to the computational procedure used 
to combine sea spectrum and RAO to arrive at a response 


soectrum. 


Development of RKAO's from Model Data 


The most comvrehensive source of destroyer model 


motions data available at the present time is the Dutch 





data that can be found in (2). It is particularly 
unique in that model response at various headings to 
the waves was measured and this is, of course, a 
necessary prerequisite if one desires to make predic- 
tions in short-crested seas and to investigate the 
general effect of heading on ship response. This Dutch 
Gata is also unique in that phase angles were measured 
and presented. ‘The primary limitation of this data is 
that the above variations were done at only one speed, 
namely 17 knots full scale. 

The vertical acceleration at the stern of a ship 
is caused primarily by pitching and somewhat less by 
heaving. Roll is the only other possible contributor 
to vertical acceleration, but its effect was ignored 
because the stern VIS towpoint is on the ship's center- 
line. The vertical acceleration also is a function of 
the phase angle between heaving anc pitching. Conse- 
quently, the data from (2) used in this investigation 
consisted of heave, pitch and the heave-to-wave and 
pitch-to-wave phase angles. Samples of the model data 
are shown in Figures 1 and 2. 

The model used was 17004. It is a model of 
Her Netherlands Majesty's Ship Friesland. Its hull 


parameters are as shown in Table I. 
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TABLE I 


L 570 ft. Waterplane Coefficient 0.805 
Beam 38.6 ft. Prismatic Coefficient 0.688 
Draft 13.2 ft. Displacement 5010 tons 
Block Coefficient 0.566 Number of Rudders 2 


It can be seen from Figures 1 and 2 that a certain 
amount of extrapolation of the model data was necessary. 
On the high end of the A/L scale, 2,>h and ¥~°2Nh/ = 
1 but the question immediately arose as to just 
where on the @/L scale it would be reasonable to assume 
that this occurs. Extensive examination of other destroyer 
model data (3), (4) showed that a @/L of 3.75 to 4.00 is 
reasonable for heave and pitch having attained their 
limiting values. Since the angle between ships course 
and wave direction had to be taken into account, Yo was 
faired to effective wave slope (based on effectiveA, 
where effective Ais the Aapparent to a ship proceeding 
obliquely through the waves) and Zo was faired toh at 
whatever actualAresulted (about 3.00 in the cases shown 
in Figures 2 and 3) when effective A/I73.75 to 4.00. 

As for phase angle extrapolation, the statements 
in (2) are generally supported by vibration theory for 
a single degree of freedom system. The limiting values 


to which by and o were extrapolated (using the same 


effective 7 limits as described above) were taken from 

















the following: 


4 ranges from (-180° to 0°) in head, 
- pow and beam seas as effective A/L 
increases from zero to about 4.00. 


5 ranges from (0° to 0°) in quartering 
z  geas over the effective A/L range of 
zero to 4.00. 
é was zero throughout the A/L range for 
2 following seas because of very small 
encounter frequency. 
ranges from (+90° to = 90°) in head and 
bow seas as effective A/L increases from 
zero to about 4,00. 
ca was 90° throughout the A/L range in 
beam and following seas. In the beam 
seas,A >> ship beam and thus there 
is effectively zero pitch phase lag 


and in the following seas caseW, is 
very small, 


v4 ranges from (90° to 90°) in quartering 
SEAS 


Before proceeding further, it is in order to pre-e 
sent the basic conventions and relationships used so as 
to prevent confusion. There are many different systems 
in use so the best policy obviously was to set up a 
system and adhere to it. I have followed that policy 
based on these conventims and relationships:* 


(a) If the wave motion amidships is 
described by 


h =h cos (, t) 
the resulting pitch motion is described by 
Ys = ‘te cos (We t +£,,) 
#In the notation used here a bar (~~) above a symbol im- 


plies an amplitude as distinguished from an instantaneous 
value. 











(b) 


(c) 


(a) 


(e ) 


(f) 


(g) 





and the heave motion by 


Z = Z, cose t + 22) 


heave is taken as positive upwards 
and 

positive pitch angle is teken as bow 

downwards. 


it can be seen from the equations in 
(a) that negative phase angles imply 
motion lagging wave and positive 

phase angles imply motion leading wave 


The phase angle between heave and pitch 
is defined by the equation 


an é- ey 


Negative values of 2, imply heave ' 
lagging pitch end positive values 
of 3 imply heave leading pitch. 


From basic wave theory we inmow that 
maximum wave s lope 


Uinax at 
It can be shown that encounter 
frequency 


! 


le = ae Ys cos, % 
g 


where “= J21g/, according to basic 
wave theory. 


The Dutch models were tested in waves 
of proportions ; 


2h = 
a 


Equipped with (a) through (g) above, it is 


possible to determine the amplitude of vertical dis- 


placement at any 


“See Figure 3 for angular nomenclature. 
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point x along the ship's length, 4,, from 

By” = Bo” + x? i? _ 28% Yocom & ey (5) 

where x = distance along snip's length from amidships to the 
point under consideration, (+) forward (-) aft. In this 
thesis x = -181 feet for the 370 foot ship length and x 
changes ina geometrically similiar manner for the other 
ship lengths investigated. The x distance was chosen 

to be a reasonable location of the VDS towpoint. 

The amplitude of vertical acceleration at any point 

x along the shipts length, Ax, can be computed from 
AX = Wef 2x 
The ordinate of the RAO is 
Beir] 
1 Ay 
and tne abscissa is ln wW,. 

RAO's for ff = 180°, 135°, 90°, 45°and 0° were 
computed directly from the model cata curves, extrapolated 
as necessary to cover a ln W range of +0.2 to -1.6. The 
reason for the limits of this range will become evident 
later. 

Then, in order to have RAO's at 22.5° increments, 
the original RAO's were first plotted on an W base to 
align the peaks and then crossplotted at constant values 
of We. The resulting RAO's for gt = 157 1/2°, 112 1/2°, 

67 a /2° Sandee ly were then converted to a ln W base 


and plotted with the original ones, resulting in a 


186, 
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complete selection of RAO's at 22 1/2° increments as 
shown in Figure 4. 

The RAO's thus derived are for the full scale 
length of the destroyer hull considered, namely 570 feet 
(this was the L used in the A/L scale of the original 
model data). 

To arrive at RAO's for geometrically similar 
hulls of other lengths is a simple matter when they are 
plotted on a ln W base (5). £As stated previously, 

w = Vene/a 
At ”~/L = a constant value, 


W, = kVl/by Subscript 1 indicates 
original value 


Wb = k Y1/Le Subscript 2 indicates 
changed value 


and 


> =| Ie/ty 


Taking the natural base logarithm of both sides 
of this equation, 1n Wy - In Wp = 1/2(1n Ig - In Lj) ana 
therefore, ln Wo - ln W, = S in @ = -1/2(1n le - 1n Lj). 


-1/2 in (78 


It follows that: dinw 


This last equation shows that RAO's for other 
ship lengths can be obtained by translating the original 


RAO's along the abscissa. 


Le 








EPHRATA PO 


a 
~Gies 
’ 3 
. 
e : 








For example, let din W = 0.1 


0.1 = -1/2 In Ig 
1) 
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and Loe = .819 (370) = 302! 


Similarly, if f1n w 


tl 


-O.1 
: = 1.222 and Lo = 452! 


aT} 


Thus, if the RAO's are shifted to the right, keeping the 
ln ® scale fixed in space, ship lengtn increases and 
conversely a leftward shift causes a decrease in ship 
length, Shifting increments of 0.2, 0.1, -0.1; and 
-0.2 were chosen, resulting in RAO's for snipv lengths 
me c29 Ttw 502 Ttege 292 Ft. and Sde2 ft. 

Figure 5 shows the 157 1/2° RAO shifted as 
described by the above procedure. All the others were 


shifted in exactly the same manner. 


Selection of Sea Spectra 


Through the principle of superposition, the 
large number of component waves comprising an irregular 
long-crested seaway may be represented in a spectral 
form having an ordinate representing energy and an 
abscissa representing wavelength or frequency. This is 
a well-established concept and is described rather well 


in (6). (lhe term long-crested is used by oceanographers 


i4 




















to describe the situation where all the waves in a system 
are coming toward the observer from the same direction 
such that he sees a welledefined pattern of long, parallel 
crests). 

The log-slope form of sea spectrum representa- 
tion developed by Lewis (5) was used in this present 
work. For my purpose it has one large advantage over 
other representations which have, for example, (wave 
amplitude )* as ordinate and encounter frequency or 
lag number as abscissa, namely: 

RAO's for geometrically similar hulls 
of different lengths can be easily 
developed by the simple expedient of 
translating them along the linw axis. 

The ordinate of the log-slope sea spectrum is 


2 
E (In & and the abscissa is inw, 
A/en 


Since a statistical sampling approach was de- 
Sired for this analysis, the Pierson-Moskowitz collec- 
tion of sea spectra (7), (8) was chosen because it is 
the largest, most recent collection of actually measured 
sea spectra. Pierson and Moskowitz use "corr .ft°"% as 
their spectral ordinate and lag number, H, for their 
abscissa. See Figure SA for an illustration of a typical 


Pierson-Moskowitz sea spectrum. The pertinent conversion 





#Corr.ft” means "corrected" ft©-Pierson and Moskowitz cor- 
rected some spectral ordinates to smooth out occasional 


unreasonable values due mainly to vagaries in the measured 
waveheight data. 
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relationshins are as follows: 


H = 180 x frequency =180~” 
eT 


‘ va 5 2 
r(inw) es [rc )| = Ww? x 180 x (corr. ft”) 
beau i er oe = 


A short computer program, written in Fortran II 


I 


source language, and shown in Figure 6, was used to cone 
vert the.Pierson-Moskowitz Spectra to the log-slope form 
used in this thesis. A sample of tne print-out is 

shown in Figure 7. 

The next matter of interest is that of shorte- 
crestedness. This is not an easy matter to dispose of, 
Simply because a snortecrested sea may arise from 
different causes. As its name implies, a shortecrested 
sea is one where the waves are coming toward the obser- 
ver from various directions, thus producing a shorte 
crested, choppy effect. The difficulty arises from the 
fact tnat a shorte-crested sea can be caused by one 
storm system, two storm systems or any number of storm 
systems. If we placé ourselves on a ship lying to in 
the ocean, maintaining a steady true compass heading, 
it is perfectly possible to conceive of a storm system 
a few miles off the bow sending wind waves toward us 
and another storm system a long distance off the star- 
board quarter sending swell waves toward us. But there 
could just as well have been just one storm system or 


more than two and they could have been off in many com- 
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binations of distance and direction, so it is only 
possible to express what is happening at a certain 
place in the ocean at a certain time by a forecasting 
or hindcasting technique. Probably the most typical 
situation is a combination of wind waves caused by a 
local storm and swell coming from a distant storm or 
storms. This is a matter of great current interest to 
oceanographers and it is my hope that in another five 
to ten years they will have accumulated a large enough 
body of directional sea spectra so as to make a reae- 
sonably sized statistical sample possible. 

However, for the moment, we must content our- 
selves with something less in our prediction methods, 
namely, either the use of long-crested seas or "uni- 
storm" short-crested seas. It is the latter that I 
have chosen. This is a situation where one storm dis- 
turbance, because of its movement anc irregular winds 
causes some spread in the wind direction and thus leads 
to a certain tyve of short-crestedness in the seas. 
These seas are characterized by a large amount of energy 
coming from a certain ee tion and progressively lesser 
amounts of energy as we look away from that direction 
until there is essentially zero energy coming from about 


70 degrees on either side of the original direction. 


aL 








This directional distribution of seaway 
energy may be characterized by multiplying a long- 
crested spectrum by a "spreading function." The one 
I have chosen is based on the work of Pierson, Neumann 
and James in (6) and is 


2 cos“ w 
7 


The result of this directional spreading of 
seaway energy can be seen at the top of Figure 9. I 
have chosen this technique because it is computat lonally 
practical and is at the very least, one step closer to 
reality than assuming long-crested seas. It is 
pertinent to note that my experience while working 
at the David Taylor Model Basin showed that short- 
crested predictions by this method are consistently 8- 
9% lower than long-crested predictions of ship motion; 
at least for the head seas case, 

With the foregoing background in mind it is in 
order now to proceed with a discussion of considerations 
that led to a final choice of sea spectra to represent 
the sea severity range of interest in these predictions. 

It seems reasonable to state that as far as 
Ship response is concerned, there are two characteris-~ 
tics of a seaway that are particularly important, namely 
the total energy, which is proportional to the area 


under the spectrum representing that seaway, and the 
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distribution of energy among the various component 
waves, an indication of which is given by the shape 
of the sea spectrum. One need only leaf through (7) 
to see that nature provides us with sea spectra of 
many different sizes and shapes. 

Now, if we wish to investigate the effect of 
sea severity on ship response, it is necessary to 
establish a primary weather severity parameter. But 
am doing so, cautionmmust be exercisec in orcer taat 
other parameters affecting ship-response are not 
ignored. 

I have chosen significant wave height (Hy /3) 
as the primary sea severity parameter because it is 
proportional to sea spectrum area and thus total sea- 
way energy. (Hi /s is, by definition, the average 
height of the one-third highest waves in a seaway 
mecord, (6). But, in order not to neglect the effect 
of sea spectrum shape on ship response, a family of 
sea spectra was needed where at eacn increment of the 
primary sea severity parameter, Hy /3 5 there would be 
several spectra of close to the same fy /3. but having 
typical vatiations in shape among them. 

The Moskowitz family of mean synoptic sea 
spectra found in (8) is the best family presently 


available most closely satisfying the above requirements. 


CO 











The family consists of five mean sea spectra, each 
spectrum being the mean of 8 to 14 component spectra 
as shown in Table II. See Figure 8 for a plot of these 


five mean sea spectra. 


TABLE II 

Moskowitz Synoptic” Number of Component Average of Component 
Mean Spectrum Name Spectra Spectra Hi/z3 in.ft. 

2OKT 12 744 

25 8 15.1 

30 12 20.4 

o5 8 25 8 

40 14 OL ek 


In each subset, the Hy /z of the component spectra 
vary in most cases by no more than a few percent, thus 
assuring that the primary differences are those of 
spectral shape. 

By comparing the ship response caused by the 
Hy /s = 7.4 ft. mean spectrum with the response caused by; 
say the Hy /s = 20.4 ft. mean spectrum, one can determine 
the effect of total seaway energy on response. On the 
other hand, by comparing the ship responses due to each 
of the 12 component spectra in the Hy /s = 7.4 ft. sub- 
set, for example, one is able to evaluate the effect of 


spectral shape on ship response i.e., determine the 





“Moskowitz classifies his spectra on the basis of wind 
speed in knots. 
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variance in response for a particular Hy /3- To my 
knowledge, it has not been attempted in previous work 
to separate these two effects and it is my hope that 
such an an&lysis will shed some light in our quest to 
learn more about the nature of ship response to a@ sea= 
WAY. 

The Moskowitz Synoptic Family of Mean Sea ~ 
Spectra converted to the log-slope form is shown in 
Figure 8. It is noted that Moskowitz required that 
the sea records chosen for his family satisfy certain 
synoptic criteria as regards non-transient wind condi- 
tions, duration, fetch, etc, all of which tend to insure 
that the seas represented are fully developed or close 
to it. He has also eliminated swell from all records 
used. So one must keep in mind that the ship response 
predictions correspond to the nearly fully developed 
wind sea with no swell situation which, of course, is 
not a precise duplication of nature. 

When there is available a larger source of sea 
spectrum data, 1 believe the soundest approach would 
be to pick out all those records having very close to the 
same Hi /3 without requiring that any synoptic criteria 
be satisfied. Thus, if the sample were large enough, 
one would be reasonably assured that a representative 


variation in degree of sea development would be 
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included. Bach sample of records having the same Hy/3s 
then, would have the variation in spectral shave pre- 


viously referred to, and then one could proceed as in 















this analysis. 


Computation of Response Data 


Once the prerequisite information, namely sea 
spectra and RAO's are in hand, one may proceed to cal- 
culate ship response. 

A verbal description of the computations will 
be given to accompany the graphic illustration given in 
Figure 9, First, the long-crested sea spectrum is 
"spread" or directionalized using the 2 cos® Mw spreading 
function previously described. This operation results 
in the short-crested sea spectra shown by dotted lines. 

Next, the directional RAO'ts are plotted on the 
same lInw abscissa and a decision is made as to which RAO's 
are going to be multiplied by which sea spectrum component 
at each increment of Inw. Figures 10 and 10a will aid 
in this "matehing" process. Say, for example, we want 
to compute acceleration response for the head seas case 
(47 = 180°) as is shown in Figure 9. We align the heavy 
wind arrow on Figure 10 with the Head Seas marker on 
Figure 10a and see that the 180° RAO must be multiplied 


by the 0° sea spectrum component, the 157 1/2/ 202 1/2° 
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RAO's (which are the same, of course, by symmetry) 


muss be multiplied by the 22 1/2° sea spectrum com- 























ponent and so on for as many of tne RAO components as 
are affected by the various sea spectrum components. 

This is done at each increment of lInw and re- 
sults in the nine response spectrum components (only 
five curves show because of symmetry). 

Now if one imagines a third dimension to the 
Response Spectrum Comvonent curves whose axis is normal 
to the plane of the paper anc whose coordinate is “, 
it is possible to visualize the integration with respect 
to “4 of the points from the Response Spectrum Component 
curves on the crossplot for each Inwincrement to arrive 
ata point on the Integrated Response Spectrum for that 
value of Inw. After the Integrated Response Spectrum 
is thus comruted, one may integrate over lnw to amive 

2 
—t the area uncer it which is the value of (SE) for the 
case considered, namely head seas, ship length = 370 ft. 
(VE is a ship response parameter that is discussed in de- 
Meet in Appendix A. Briefly, the YE parameter reflects 
certain statistical proverties of a short-term sample 
record of response Guring invariant conditions of weather, 
Speed and heading.) The angular integration and then the 
inw integration were done using Simpson's Rule in the 


computational scheme employed. Nineteen Inwincrements 


ol 





of O.1 each were used along with the 22 lye angular 
increments as shown. The risht and left hand extremities 
of the Inwscale, namely, -1.6 and +0.2, respectively, 
were determined by the termination of most typical sea 
spectra on the right hand end and the termination of the 
RAQ'ts on the left hana end. 

The process described above can then be done for 
the eipht other headings between head and following seas 
so that each sea spectrum and each set of RAO's for a 
given ship len.th results in nine JE values, one for 
each heading from head to following seas at 22 1/2° 
increments. 


This was cone for the five Moskowitz Synoptic 


ean spectra for each of the fiwe mip lengths and at 


nine headings, resulting in the 225 \IE values found in 
Ficure ll. 

Since each one of these VE values is the result 
of the process of computation just described, where sea 
spectrum components are multiplied by RAO components to 
arrive at the Response Spectrum Components which are in- 
tegrated to arrive at the Integrated Response Spectrum 

hich in turn is integrated to arrive at the VE value, 


it can readily be seen that a very large number of 


numerical calculations are involved. However, since 
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these calculations are of a repetitive nature, including 
several cyclic permutations, they are especially amenable 
to calculation by electronic computer, and this is how 
they were carried out. 

The program shown in Figure 12, based on a 
similar one in use at Webb Institute, was written by the 
author in the Fortran II source language while on temporary 
duty at the David Taylor Model Basin and the calculations 
were carried out on the UNIVAC LARC facility there. A 
sample of the print-out cm be seen in Figure 13. It 
might interest the reader to know that the computer took 
less than three minutes to compute the 225 VE values and 
having done similar calculations with a desk computer, 


LI estimate it would take about a month of 8-hour days 


to calculate the same number of E values manually. 


Thus it is seen how the electronic computer cm greatly 
enhance the scope of an investigation like this one. 

The tabulation of VE values in Figure 11 shows 
the effect of seaway total energy on acceleration rese= 
ponse, but it does not show the effect of sea spectrum 
shape . 

As previously described, the effect of sea spectrum 
shape on ship response was to be evaluated by computing 


\/E values for each of the component sea spectra in the 
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five Moskowitz subsets. Since there are a total of 54 
component spectra in the Moskowitz Synoptic Family of five 
mean epectra, converting each of these to the log-slope 
form, plotting them, tabulating and card-punching them 
as input to the ship response comouter program would 
have been too formidable a task to undertake in this 
thesis. Instead, the above was done only for the 12 
component spectra of the Hy /3 = 7,4 ft. mean spectrum 
and the 14 comynonent spectra of the Hy /s = 52.2 ft. mean 
spectrum so as to bracket the range of tne primary sea 
severity parameter Hy/3- The resnonse characteristics 
for intermediate values of Hy /s were obtained using 
straight line interpolation which turned out to be a 
culte reasonable procedure as will be discussed later. 

The VB values for all nine headings and five 
ship lengths in the 12 component sea spectra of the 


Hy /3 = 7.4 ft. subset and in the 14 component sea spectra 


of the Hy /s = 32.2 ft. subset were computed using the 


same computational scheme as was previously described and 
then for each combination of ship length and heading, the 
following important characteristics of the VE distribu- 
tion were cetermined: 


(a) the sample mean of \F was computed using 
m™ 


oS : d nat 


where n = 12 or 14 depending on which subset 
was under consideration and X = sample mean 
JB and x1 = (VE), 








(b) the variance of YE was computed usiny: 


Cr J (xs )% - i°| 
S 1 (¢ Soe n = 


again where n = le or 14 and s© = sample 
variance and s = sample deviation 


These conputations also were done on the DTIMbB 
UYIVAC LAKC using the Fortran II source program written 
for this purpose which is shown in Figure 14. This 
program is similar to the one in Figure 12 except that 
a section for computation of sample means and devia- 
tions was added onto the section that computed the VE 
values. The print-ovt of results from this program 
may be seen in Figures 15 and lda. 

At this point, one small issue was disposed of. 
For each ship length and heading, the mean WE values in 


Figure 15, for example, should agree with the VE values 


Seven in Figure 11 for the Hy /s = 7.4 ft. mean sea 


spectrum, that is to say that the mean value of VE 
obtained using the subset component spectra should be 
equal to the WE obtained Using the mean sea specurun. 

This is mainly a crosscheck on the computational scheme 
used. A comparison of the data shows a maximum difference 
in the order of 4 percent, so the procedure of using a 
mean sea spectrum to obtain mean UE response seems 
justified. This conclusion engenders confidence in the 


JE values shown in Figure ll, because the ones for the 
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three intermediate values of Hy /3 were computed using the 
mean sea spectra only. 

The sample JE deviations were linearly inter- 
polated over the Hi /3 range and are shown in Figures 
16 a-e for the various ship lengths and headings. 

To summarize, it is seen that at this point, 
enough information is available so that the effect on 
acceleration response of both sea spectrum total energy 
and sea spectrum shape can be examined. 

The effect of a variation in sea spectrum total 
energy manifests itself in the five groups of sample 
mean \E values in Figure 11. The effect of a variation 
in sea spectrum shape manifests itself in the sample 
VE deviations shown in Figures 16 a-e; these deviations 
from the mean exist solely because of this shape effect 
on ship response. 

At this point, the WE values in Figure 11 could 
be variously plotted to show the trends of VE stern 
acceleration with heading, length and sea severity, but 
before this is done, it is desirable to treat these VE 
values which are characteristic of short-term response, 
as a statistical sample and use them to achieve a reae 
sonable estimate of what SE values might be expected within 
the total universe or population of all possible destroyer 
stern acceleration \E values during a destroyer's lifetime 
under similar conditions of heading, length and sea severity. 


The next section will be devoted to this estimation. 
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STATISTICAL ESTIMATION OF POPULATION YE VALUES 


me eee DR ee 


As discussed in the previous section, at each 
combination of heading,ship length and Hy 3, there is a 
mean value of UE (Figure 11) and a distribution of YE 
values above and below that mean, wnich distribution is 
characterized by the sample deviation found in Figures 
16 a-e for the same combination of heading, length and 
Hy /3. 

The assumption was made thet this distribution 
is Gaussian or normal. The normal probability distribu- 
tion is completely defined by a mean value and the de- 
viation from that mean. the primary defense of this 
assumption of normality is that it does not seem un- 
reasonable in the light of present knowledge. Perhaps 
vhen much more work in this area has been done, a slightly 
better distributim can be found. ‘The size of the statis-~ 
tical sample used here is too small to permit construction 
of a valid histogram. This assumption of normality has 
also been used by other investigators (12). 

lhe next and last hurdle is how to use a sample 


mean and deviation to predict the true total population 


2 








mean VE and its deviation. This is where one may make 
good use of the concept of statistical inference, i.e., 
one may infer from the characteristics of a statistical 
sample the characteristics of the total population. In 
this case the assumption was made that in the sample, 
the VE values are normally distributed with a mean, X, 
and a deviation, 8. 

Sample: N (x,8°) 
In the population, the YE values are also assumed 


normally distributed, with a meanfand a deviation, 6. 


Population: N(M,o°) 


The task now is to estimate /7 and G on the 
basis of X and s. Presently available theoretical 
statistics gives one two choices that are quite differ- 
ent in concept, namely the processes of point esti- 
mation and interval estimation. 

I will first discuss point estimation. Maximum 
likelihood estimators may be derived which give a con- 
sistent,unblased estimate of the population parameters, 
In the case under consideration, the unbiased maximum 


likelihood estimator for Ais: 
med bx =k (9) 
n 


and the unblased maximum likelihood estimator for 0 is; 


= ah a ==)” = % 
GF ie S (x ) S 


*It can be shown that this expression is equivalent to 
the one used on page40. 











This means that the best estimate one can make 


in this case, of “and 6 is to use XK and s, respectively. 
However this somewhat vague approach is not very attracte 
tive unless the sample is of a respectable size, which 

in this investigation it is not. Perhaps the largest 
drawback to the use of point estimation is that when 
Eampilé size is relatively small, one is provided with 

no specific expression of confidence that can be had in 
the estimation of the true parameters Mand 6. 

It is due to the above that 1 decided to use the 
concept of interval estimation. With a specific level 
of confidence in mind, say 90%, one may predict a para- 
meter interval or region which, withMl 90% confidence, 
will cover the true parameter value or values being 
sought. 

It has been shown (9) that the statistic 


x -/ J 
o/un 


is normally distributed with zero mean and unit variance, 
and that the statistic 


ns@ 5 
TE 
is chi-square distributed, where So and “> are the ‘rue: 
population parameter values. 

Thus, using any tonvenient set of probability 


tables, the statement may be made that ’ 








* These two statistics are stochastically independent. 
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Pr /-1,15 y X -o 7 1.25 | = 875 
Foie 
from the standardized normal distribution and that 


pr | 4.87 L ns® é 19.7 | = .900 
Jo 


from the x distribution with n-l degrees of freedom 
where n = 12 in this case. (n = 12 was used for the sake 
of uniformity even though the Hj /3 = 32.2 subset had 14 
component spectra-this will yield slightly conservative 
results ) 

The joint probability is (9) 


P “1.15 ¢ X - %0 < 1.15, 4.57 a ns® (19.7) = 
| do/Yn cor 


(,.875)(900) = .788 
These inequalities define a confidence region in the 
Oss parameter space as shown in Figure 17. 
ca There is a 78,8% chance that the crosshatched 
parameter region covers the true parameter point (Ao, 60), 
wa wherever it may be, that is to say that if 100 samples of 
2 each were chosen _ if the 100 resulting 78,8 con- 
fidence regions were plottea as above, about 79 of them 
would cover the true parameter point and the remaining 
21 would fail Lc We so. This is the quintessence of 
the concept of interval estimatim. 
. 


i= The approximately 80% confidence level was chosen 


as being reasonable and the following discussion is 
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intended to show that it is more than conservative 


enough for the purpose of these estimations. 
The maximum likelihood estimators tell us that 
P is the most likely location of (ov, So). With this 
and the approximately 80% confidence level in mind, I 
have chosen (Mu, Gu) which will give the "worst" normal 
distribution in terms of extreme values as the true 
parameter point. This is a fairly conservative choice, 
because of the likelihood that (60,/o) is somewhere 
near P. I feel this is a situation where a little 
judicious conservatism is better than being inordinately 
bold in ignorance. 
The parameters of N(/u, 6 u”) may thus be evalua- 
ted from the followtng simple relationships: 
Au 
Cu 
These equations show the amount by which the 


x + .540 8 


1.620 gs 


sample parameters are increased so as to arrive at an 


estimate of the population parameters. 
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INVESTIGATION OF TRENDS IN VE ACCELERATION 


From this point, mean YE acceleration will imply 
estimated population. mean YE acceleration (Mu) and VE 
acceleration deviation will imply estimated population 
deviation of VE acceleratim (6u). 

From the data in Figure 11 (x) and Figures 16 
aee (s), many combinations of ship heading, length and 
Hy /3 can be used to investigate trends in VB accelera- 
tion. For purposes of the present work, I have considered: 


(a) one ship length (370 ft) - all naad aay 
all 1/3 


(b) one heading {head seas) - all lengths- 
all Hy /3 
These cases show the effect of ship length and heading 
as well as Hy /g on mean \E acceleration response and its 
deviation. Other cases could, of course, be considered 
if desired, using other data from Figures 11 and 16 a-e, 
Figures 18 and 18a deal with the 370-foot ship 
length case and show how mean VE acceleration and its 
deviation vary with heading and Hy /3- Figures 19 and 19a 
show trends with ship length and Hy /3 for the head seas 
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Finally, by crossplotting Figures 18 and 18a 
at discrete values of Hi/3 (which can be thought of as 
different sea states - see Appendix B), Figure 20 was 
derived. Similarly crossplotting Figures 19 and 19a led 
to Figure 21. In the interest of clarity, only enough 
deviations are plotted on Figure 20 to demonstrate the 
fact that the deviations hehave in much the same fashion 
as do the mean VB values for any given ship length, 1.e., 
a plot of deviations results in curves similar to the 
curves of mean YE acceleration in Figure 20. 

If one accepts this conclusion of similar be- 
havior of means and deviations, the fairly mild curvae- 
ture of the mean VE curves between Hy /3 = 7.4 and 32.2 
feet in Figure 18 tends to justify the assumption of 
linearity dae in interpolating the sample deviations in 
Figures 16 a-e in the same Hy /s range e 

There are some important points to be learned 
from an examination of Figures 20 and 21. 

Figure 20 illustrates a point of great opera- 
tional significance. It is seen that if destroyers were 
designed so as to tolerate beam, quartering and following 
seas operation from a rolling anc maneuverability stand- 
point, the stern accelerations could be greatly diminished 
and thus the effectiveness of the DD/VDS system is 


greatly enhanced in any given sea state. A step in the 
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right direction can be taken in the installation of 
anti-rolling fins on destroyer-type ships, because as a 
consequence of being able to tolerate beam seas opera= 
tion, Figure 20 shows they could diminish stern accele- 
rations to between 50 and 60 percent of the head seas 
accelerations. 

As a slight digression, I see no reason why 
the same reasoning cannot be avplied to problems associa- 
ted with bow movement such as slamming, shipping of green 
water, degradation of bow sonar performance and general 
ability to maintain high speed in rough seas. 

Figure 21 illustrates the fact that ship length 
has a greater percentage effect on stern acceleration in 
small waves than in large ones. Consequently, increased 
Ship length pays significant dividends during the majority 
of the time because the small waves are much more prevalent 
than the large ones, but pays lesser dividends when con- 


sidering relative performance in severe waves. 











UTILIZATION OF xf ACCELERATION DATA™ 


el es OEE ee ee Eel 


A crucial question is how to most effectively 
make use of data such as that plotted in Figures 20 and 
el. By the means discussed in Appendix A, it is possible 
to predict from the JE value the probability of occurrence 
of any particular value of stern acceleration. 

| However, it is highly desirable to simultaneously 
consider the normal distribution of VE values (and thus 
the Rayleigh parameter R) and the Rayleigh distribution 
of acceleration single amplitudes. The theory for this 
has been developed by Nordenstrom in (12). He assumed 
that the parameter R in the short term Rayleigh dis- 
tribution is normally distributed and proceeds to obtain 
overall distributions by numerical integration of the 
convolution of the Rayleigh and normal distributionse A& 
sample of the results of this numerical integration can 
be found in Figure 22 with suitable nomenclature changes 
made to adapt it to the present work. One can enter this 


diagram with 6u//u and for different values of N (number 





"a detailed discussion of the VE parameter is given in 
Appendix A, 
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of cycles) obtain a factor, which when multiplied by Mu 
allows one to derive curves of Ax/g such as those shown 
in Figure 23 for the head seas and 370 ft ship length 
CASE. 

It must be pointed out here that although these 
curves in Figure 23 could be termed shorteterm distribu. 
tion curves in accordance with Jasperts definition (11), 
I will not do so, because it is absurd to consider that 
millions of cycles could take place during a short term 
of say, twenty minutes. However, I do not choose to call 
them long-term distribution curves, either, because this 
implies a natural vaiation of ship speed and heading, 
which is not the case here. These curves will be con- 
sidered "long-term" only in the sense that each curve 
is the sum of many short-term distributions, keeping in 
mind the assumptions of invariant ship speed and heading 
that I have been forced to make, the former because of a 
paucity of model test data and the latter because of a 
lack of information available as to probability distribu- 
tion of heading to the sea. This matter is discussed 
further in the next section. 

A partial answer to this dilemma is to consider 
the probability of occurrence of the five sea states 


shown. It is then possible to weight these five curves 
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by the probability of occurrence of each sea state and 
arrive at a single curve that is at least one step 
closer to being a true long-term distribution curve. 
The following sea state distribution from (14) 
for the North Atlantic (40-50° North Latitude) was used 


as an example: 


Hj /3 in feet Probability of Occurrence 
Calm 1.3% 
less than 4.92 Sea? 
4,92 - 9.84 44.1 
9,84 2 14.77 14.9 
14.77 _ 19.69 509 
19.69 = 26.20 2 
greater than 26.20 Lee 
00. 


Plotting this information on probability paper 
gives one a continuous distribution of Hy /3. By breaking 
this continuous distribution up into discrete increments 
centered on the Hy /3 values being used for prediction 
purposes (6,11,17,25 and 37 feet), it is reasonable, but 
not mathematically rigorous to state that: 


Ha /5=6 ft is representative of what occurs 70.0% of the time 
tt 


ot : 22.0 . 

17 2 5.5 7 

! 25 "! 2 6) | 

it a7 hy "5 iy 
100.0% 


These statements are not mathematically precise 
because the ship response due to a discrete Hy /s is said 
to be equivalent to the ship response resulting froma 


range of Hy /s values above and below that discrete value, 
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and this is obviously affected by the size of the Hy /s 
increment between discrete values of Hy /3+ The larger 
the increment, the larger would be the deviation in VE 
values; larger than the deviations caused merely by 
variation in sea spectral shape. 


However, in the analysis undertaken here, the 


ee ee ———UlUl”UlUCOCOOOOES 


degree of approximation seems reasonable with general 













regard am to accuracy of the predictims. 
The five curves in Figure 23 were weighted at 
. various levels of acceleration by the above probabilities 
and the result was the curve shown in Figure 24 which I 
choose to call a quasi~long term distribution curve for 
a samt veGlat 1eRceh ship under the assumed conditions 
of ship speed and headtne, When statistical information 
is available regarding speed and heading distributions, 
inclusion of this information in the analysis will make 
@ curve such as in Figure 24 completely general, and 
then it could honestly beccalled a long-term distribu- 
tion curvese* mp 7 
Fortunately, however, a quasiwlong term distri~ 
bution curve such as that in Figure 24 is presently 
useful for systems anduysis and design because the oon- 
ditions of apeed and heading used in this work are such 
as to yield conservative results. In other words, more 
work and increased knowledge will show us how much 
Lower our predictims can safely be. Here again this 


_ is given more disousgion in the next section. 
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. A DEMONSTRATION OF THE SYSTEMS VIEWPOINT 


With predictions available such as those 
developed in the preceding section, it is possible to 


make an evaluation of system effectiveness, keeping in 





mind the assumptions as to ship speed and heading that 
" were made. The significance of these assumptions will 
be discussed further later. 
; In the case under study, the system is the 


DD/VDS system. Destroyer vertical stern acceleration 












is no doubt the most important single parameter affecting 
Operation of the DD/VDS system, given an otherwise oper- 
able system. The structural ability of the cable to 
remain intact in rough weather is a go/no=go situation 
Basin adlott on vertical stern acceleration. Vertical 

_ stern acceleration also bears on understanding and 
solving the VDS control problem and on the design of 

the towing device. 


> - A systems analysis has been done herein using 
vertical stern acceleration as the significant parameter- 


: other parameters could be included in the same manner, 


. “ae 
. «af such wie deemed necessary. 
‘ 7 


ve 
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There are two basic areas of interest in this 











analysis of DD/VDS operatimal effectiveness based on 
vertical stern accelerations; they can best be posed as 
questions: 

(a) in the case of an extant system, where a 
design assumption as to vertical stern 
acceleration has been made, what should 
the upper operational limit be, with 
respect to weather, in order to give 
reasonable assurance that the cable will 
not be snapped and that the VDS might be 
controllable, that is to say, in order 
that the design assumption of vertical 


stern acceleration will not be exceeded? 


(bo) in the case of a system under original 
feasibility or design study, what value 
of vertical stern acceleration should be 
used as a design requirement, so as to 
insure that the system will be operae- 
tionally effective for an acceptable per- 


centage of the time on a long-term basis? 


For the extant system, the curves in Figure 23 are 


the most useful. Say, for example, that Ax = 1.5 g's was 
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used as a design assumption and therefore the design was 
based on this figure. Let us say, further, that it is 
desired that the probability of exceeding 1.5 g be .001% 
or less. Figure 23 tells us that in seas characterized 
by H1/3 = 25 feet, the probability of exceeding 1.5 g is 
slightly less than .001%,which is satisfactory. However, 
in the case of seas characterized by Hy /s = 37 feet, the 
probability of exceeding 1.5 g is about .095% ,wnich is 
too great. Therefore, we would specify that this system 
avoid operating in seas having Hy1/3 greater than about 
25 feet. Since seas of this severity might be expected 
to occur only two or three percent of the time, this 
system could be expected to be operationally useful 
97-98 percent of the time, considering the restriction 
imposed and the probability of failure that was 
specified as being acceptable. The numbers used in this 
example are for illustration only - if a smaller prob- 
ability of failure is desired, for example, a lower 
upper sea severity limit would have to be imposed, thus 
reducing the total percentage of the time that the system 
could be operated. This is obviously a situation involving 
various trade-offs and one must pick the combination best 
suited to his needs. 

In the case of a new design, it can be seen that 


Figure 24 is the most useful, if it is desired, as it 
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certainly seems it would be, to avoid any sea severity 
upper limitation on the operation of the system. Once 
an assumption is made as regards likely areas of opera-= 
tion, a curve such as that shown in Figure 24 can be 
used to arrive at a reasonable value of stern accelera- 
tion for use as a basisfor design. 

For example, let us consider a DD/VDS system 
that is at sea for 40 percent of a ten-year life span 
(the ship will last longer, but most likely the weapons 
system will be changed), and assume that the VDS is 
streamed for half the time tite ship is at sea. One 
could expect about 9 million cycles of stern accelerae- 
tion during that period of time, based on ™® overall 
period of stern motion of about 7 seconds. Entering 
Figure 24 with N = 9 million, one reads off about 1.8 g's 
as the level of stern acceleration above which we can 
Bzpect only one of the 9 million cycles to occur. The 
corresponding probability of failure if we use 1.8 g's 
as a basis for design is about .00002 percent. 

It must be emphasized again at this point, as in 
the previous section, that curves similar to those being 
discussed are only as good as the assumptions made in 
deriving them. In the specific case under consideration 
in the present work, the curves in Figures 23 and 24 are 


based on a single ship speed of 17 knots and, of coarse, 











the system will not spend its entire lifetime oprating 
at 17 knots, although some weight may be given to the 
fact that 17 knots could be considered a typical speed 
for operation of the DD/VDS system. Another fact that 
must be kept in mind regarding these particulor figures 
is that they are for head seas operation and, here again, 
the system will not spend its lifetime operating in head 
seas. As a matter of fact, as discussed previously, it 
would considerably lower the levels of stern accelera- 
tion if we were to design our destroyers with anti- 
rolling devices so as to be able to avoid head seas 
Operation. Therefore, the curves given here will yield 
conservative results in any case. 
My primary purpose in this section is to present 

a method for systems analysis that can be quite useful 
in evaluation of system effectiveness when the input to 
the analysis is reliable in magnitude and sufficient in 
scope. At present, I consider the available input in- 
formation marginal in some respects, namely: 

(a) probability distributions of ship 
heading to the sea ge not presently 
available and should therefore be 
given some attention. It seems evie 
dent to me on the basis of four years 
at sea on board destroyers that as 


weather becomes more severe, the 
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tendency to steer so that seas are on 
the bow increases because of a desire 
for minimum motions and maximum 
maneuverability. It is due to the 
foregoing that I assumed head seas as 
being slightly more conservative than 
an assumption of bow seas. However, 
the more widespread use of such de- 
vices as antierolling fins will no 


doubt change the situation. 


(b) it would be highly desirable to include 
the effect of speed in these mwedictions 
even though it probably has a second order 
effect compared with ship length and sea 
severity. At present there is no model 
test series available that covers all 
headings and speeds. Perhaps analytical 
derivation of the RAO's, such as is now 
being investigated at Webb, would solve 
this dilemma, by replacing the model 
testing process with an analytical pro- 
cedure than can be done with a computer. 

One further point that must be discussed is that 


of attempting a complete systems analysis. There are 
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many factors that affect a system as complex as a ship 
and its equipment, and I have examined only one sub- 
system of the larger overall system. The reliability 
of the machinery in its task of motivating the ship 
during VDS operations at the desired speed certainly 
affects the effectiveness of the DD/VUS system, but the 
inclusion of factors such as this is left to the opera- 
tims analyst charged with maintaining the macroscopic 
view, more commonly referred to as the "big picture.” 
It is my hope that the information developed in this 
thesis regarding destroyer stern acceleration might aid 


him in his task, 
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CORRELATION OF PREDICTIONS WITH FULL-SCALE MEASUREMENTS 


As a result of interest in VDS systems, the U.S. 
Navy Underwater Sound Laboratory in New London, Connectie- 
cut, has recently embarked on a program of collecting a 
statistical sampling of vertical stern accelerations on 
destroyers at sea (17). For this purpose, Giannini 
Controls Corporation Model 2432 statistical accelero= 
meters have been installed on various destroyers at the 
stern, as well as other locations. 

Up to the present time, these accelerometers have 
been installed so as to detect instantaneous values of 
positive acceleration that exceed the following four 
thresholds; 1.20g, 1.55g, 1.50¢ and 1.65g, and record 
on the four counters, the number of exceedences above 
each of the thresholds, 

it must be pointed out here that since these 
accelerometers are designed to be used horizontally as 
well as vertically, the zero reference for measuring 
vertical accelerations is one g rather than zero g, the 
latter being used throughout this thesis. As a conse= 
quence, the above thresholds become .2c0Og, 55g, .o0g 


and .65g in the zero g reference system for vertical 


PL. 











accelerations. No attempt will be made here to debate 
the pros and cons of the two systems - suffice it to say 
that there are arguments for both. 

As of the present, there has been collected a 
small amount of this statistical data (and to my know- 
ledge it is the gnly such data available) and it is of 
considerable interest to compare it with the predictions 
made in this thesis. 

The data received from the Underwater Sound 


Laboratory is as follows; 


USS MOALB (DD693) 


Number of accelerations above .20g = 708 
.O5g = 259 
00g = 44 
65g = 4 


Total time of test = 226 hours 

Estimated number of cycles for 6 second 
period of motion = 135,600 cycles 

Weather 1/6 - 6 foot wave heights observed 

Headings and speeds = various 


The probability of exceeding varias levels of 


acceleration, therefore, would be as follows: 


Probability of exceeding .20g = 708/135,600 = .530% 
.o5¢ = 2310/155,,6008= .176 
50g = 44/135,600 = .0325 
.65¢g = 4/135,600 = .00294 
USS PUTNAM (DD _ 757) 
Number of accelerations over .18g = 6009 
06g = L529 
46g = 149 
672 = LY? 
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Total time of test 612 hours 

Estimated number of cycles for 6 second 
period of motim = 367,000 cycles 

Weather 1-10 foot waveheights observed 

Headings and speeds = various 


Probability of exceeding .18g = 1.64% 
08g = 45 
46g = .041 
672 = «0046 , 


USS H. J. ELLISON (DD 864) 


Number of accelerations above .19g = leeg 
04g = 156 
48g = 18 
66g = Z 


Total time of test 629 hours 

Estimated number of cycles for 6 second 
period of motion = 577,000 cycles 

Weather O-l2 foot wave heights observed 

Headings and speeds = various 


Probability of exceeding .19¢ =.34% 


04g =.045 
-48e =.0049 
66g =.0005 


The Moale and Putnam are destroyers of the Sumner 
Class whose length on the waterline is 369 feet. The 
Bllison is a Gearing Class Destroyer wnose length on the 
waterline is 385 feet. 

The probabilities of exceeding the four levels 
of acceleration for each ship were plotted in Figure 23 (Pg.70) 
to compare them with predicted values. 

Since the predicted curves are for head seas 


and 17 knots ship speed, one would expect them to be 








high relative to the observed data which is for various 
headings and speeds. Therefore, since the observed data 
was taken in seas in the neighborhood of O-1l1 feet in 
height, one would expect the observed data to plot some- 
where below the Hy /s = ll feet curve. As can be seen in 
Figure 235, this is exactly where the observed points plot. 
The dashed line shows the trend for the two 369 foot 
ships and the dot-dash line shows the trend for the 583 
foot ship. The longer ship exhibits a tendency toward 
lewer stern accelerations which agrees with the trends 
of stern acceleration with ship length previously dis- 
cussed. 

The fact that the lines drawn through the ob- 
served data have a larger slope than the individual Hi/s 
Curves seems to reflect the fact that the observed data 
included the effect of a range of weather conditions. 
The curve of predicted values in Figure 24 has a larger 
slope and it likewise includes the effect of a range of 
weather conditions. 

This section is intended in no way to be a come 
prehensive correlation of observed data with predicted 
data. Indeed, much more observed data is needed so one 
might cull it for information as to accelerations for a 
particular heading and speed. However, the information 


presented here tends to increase one's confidence in the 
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quantitative and qualitative aspects of the stern 


accelerations predicted by the method presented in this 


thesis. 
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CONCLUSIONS 


1. The method presented in this thesis pro- 
vides the system designer with a powerful technique -— 
where ship response to the sea is a pertinent factor—- 
for analyzing on paper the worth of a proposed system while 
still in the conceptual stage.’ | 

2. The available full-scale stern acceleration 
data tends to confirm the qualitative and quantitative 
aspects of the predictions made by the method presented 
herein. 

5. From the standpoint of diminishing stern 
acceleration, ship heading to the seas is the most impor- 
tant single controllable factor and, therefore, destreyér- r 
type vessels with installed equipment sensitive to atefh: 

, accelerations should be designed so as to be able to 

. tolerate rough seas operation in beam, quartering, and 

. following seas. Installation of antd-rolling devices is 
} a step in the right direction by enabling the ship to 

: operate in beam seas and thereby paying a dividend of 50- 
, 60 percent diminution in stern acceleration. 

4. Ship length has a greater percentage effect 
on stern acceleration in light seas than it does in 
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heavy seas. Consequently, increased ship length pays 
significant dividends during the majority of the time 
since the light seas are much more prevalent than are the 
heavy ones. The effect of ship length is not so marked 
when considering relative performance on a percentage 
basis in severe seas, 

S. Sea spectrum shape has a significant effect 
on stern acceleration response and should therefore be 
accounted for in the method used to predict that response. 
The method used here seems to be effective in that regard. 

6. As demonstrated in this thesis, the system 
@pproach seems particularly useful in certain aspects of 
ship design, since a ship-related system may most logically 
be cmsidered in terms of its intended mission, for example, 
the location of enemy svibmarines,. 

7. Ina problem such as the one studied here, 
the use of a high speed electronic computer virtually 
frees the designer from worries regarding time required 
for calculations and allows him to spend full time ine 
suring the soundness of the analytical processes being 
used. 

8. The generality of the analysis developed 
here insures that a similar method can be utilized to 
investigate other system characteristics related to sea= 
keeping, such as shipping water over the bow, sonar de- 


gradation and slamming. 
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9. The analysis used in this thesis shows that 
the behavior and effectiveness of a ship-related system 
is ideally suited for expression in probabilistic terms. 
However, both ship designers and operators will need to 
develop a whole new framework of thought and reference 
as to acceptable probability levels of system failure 
and/or loss. Indeed it puts a designer face to face 
with his responsibilities when he has to express loss of 
resources, human and otherwise, in mathematical terms 


that are uncomfortably precise. 
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RECOMMENDATIONS FOR FUTURE WORK 


l. Investigate means that would allow destroyer- 
type vessels to tolerate operation in rough quartering 
and following seas. The obstacles to be overcome include 
considerations of stability and maneuverability. 

2. Continue work presently being done in the 
development of theory for analytical derivation of res- 
ponse amplitude operators, including experimernal verifi- 
cation. 

Oe Investigate probability distributions of ship 
heading to the sea under various conditions. 

4. When data becomes available, undertake a 
comprehensive correlation between full scale stern 
acceleration measurements and the predictions contained 
herein. 

5. Give increased emphasis to collecting more 
seaway measurements useful for ship design, preferably 
directional data measured with a buoy rather than a 
shipborne recorder. In this regard, effective liaison 
between ship designers and oceanographers is of the ute- 
most importance, 

6. Investigate stern vertical velocities and 


displacements under various conditims of weather, ship 
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length anc heading using a method similar to the one 
presented here. Also investigate the effect of ship 
speed on stern vertical acceleration, velocity and dis- 
placement. 

7. Investigate the application, to the VIS 
cable stress situation, of an instrument similar to the 
flexing stress monitor described in (18). This instru- 
ment could be used to compute and prescnt maximum probable 
cable stress during the ensuing 15 minutes or so of 
Operation and thus give the ship's personnel an objec- 
tive basis for making speed and/or heading changes when 


the probable cable stress becomes excessive. 
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APPENDIX A 
DISCUSSION OF THE YE PARAMETER 


A record of any wave-caused phenomenom in 


irregular seas such as stern acceleration would look 


something Like 


Srey | : 
Picasa ns a [ Le 


where xj is stern acceleration measured fren the mean 
valve at equal or random time intervals and yi is peak to 
peak stern acceleration measured wherever it occurs. 

Lt is now widely accepted in the literature” that 
the x, 's are normally distributed and that the yi's are 
Rayleigh distributed on a short-term basis, A short- 
term distribution ig applicable so long as uniform condi- 


tions of sea, ship speed and course persist. This is to 


See for example, (11) and (15). 
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be distinguished from a long-term distribution which allows 

the sea, ship course and speed to vary over a range of 

conditions. Jasper introduced these definitions in (11) 

and they are used by many investigators, including Nordenstrom (12). 


(Peak to peak) Rayleigh: f (y) = ey ‘ -y*/R 


(Sucessive - = u)* 
Ordinates ) Normal: f (x) = 1 . aly =) * 


It is seen, then, that the parameters necessary to define the 
above distributions for a particular case are R and 0 o Fortunately, 
it is possible to get these parameters from the Integrated Response 
Spectrum which, through the superposition concept, has certain 


mathematical properties in common with the irregular record shown. 


- = GREA UNDER 
SPECTRUM 


pe) 8 





* for brevity, the normal probability distribution 
function is often referred to as N(y 0). 


9b, 








It has been shown in the literaturé and is now 


widely accepted that, for the type of spectrum being used 





in this work, 















R= 86% = 45 (13) 
Thus it is seen that R and6 * are known once E is known 
from the Integrated Response Spectrum. 
It can also be demonstrated, using the concept 


of maximum likelihood estimatim, that 


R=5 ne 


Aa n 


- and since root mean square (of peak to peaks) is defined 
»* 
as 


ak to 
-" ) peaks 


_ one may relate RMS(p4p) to R and thus E as follows 


> F. 
RMS (pp) 2 | 22 = UR = 


RMS of peak to peaks is commonly referred to as RMS of 
@e 
double amplitudes. Therefore RMS of single amplitudes 


rn 


RMS/ single = RMS, /2 =UE 
amplitude PUP 


. Thus it is seen that the YE values plotted in 





the various figures can be thought of as being the root 














mean square of single amplitudes of stern acceleratim 
and such is the case throughout this thesis. 

All of the above relationships were originally 
developed using an irregular record of wave height and 
the corresponding sea spectrum, but I have discussed 
the situation in terms of stern acceleration because 
ship characteristics dependent upon waves such as pitch 
and heave have been found to behave in the same fashion 
as the waves. 

Since most often one is most interested in single 
or double amplitudes of motion, it is common practice 
to use the parameter E which defines the parameter R 
and thus the Rayleigh distribution of amplitudes to 
arrive at, for example, the average highest expected 
value of stern acceleration single amplitude in 20 cycles 
from 


( x/e)s joo = 1.87 70E 


Multipliers for other numbers of cycles N (where 
N>50) may be computed to within 3% accuracy using the 


following formula from (10); 


1/2 - 
(In N) / a: ae N) Ve where 0 = Eulers Constant = 


Pb ial 
Examples are as follows: 
=number of cycles Average of highest expected 
values of single amplitude 


50 2.12 VE 
100 2.28 UE 
200 2.43 VE 
500 2.61 VE 

1000 2.74 VE 
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For the reader's convenience, the following 


are also quoted: (6) 


Most frequent single amplitude 0.707 UE 
Average of all single ampli- 
tudes = 0.886 UB 
Average of 1/3 highest single 
amplitudes = 1.416 IE 
fore of 1/10 highest single 
amplitudes 1.800 7 











APPENDIX B 


A UNIFORM DEFINITION OF SEA STATES 

For about fifteen years, a world-wide sea state 
scale (World Meteorological Organization Code 75) (15) 
has been in general maritime use by both naval and 
merchant ships for the collection of weather data at 
sea, 

It seems logical, therefore, in the interest of 
meaningful fullescale correlation, that all our engineering 
predictions for ultimate application to the ships at sea 
Should be in terms of WMO Code 75, which is quoted below 


for the reader's convenience; 


STATE OF SEA-=WIND WAVES (WMO COLE 75) 


Code Description Height in Feet # 
0 Calme-glassy 0 
i, Calm-rippled 0-1/3 
2 Smooth-wave lets 1/3-1 2/3 
3 Slight 1 2/34 
4 Moderate 4-8 
5 Rough 8-15 
6 Very Rough 13-20 
7 High 20-30 
8 Very High 50-45 
9 Phenomenal - over 45 


"¢ is proposed herein to consider this column as being 
1/3. 
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In spite of the above, other sea state scales, 
such as the one reproduced in Figure 25, are in wide use 
by designers, presumably because of the mathematically 
convenient use of significant wave height (Hz /3) versus 
sea state designations. However, one simple and defensible 
assumption would remove the need to use anything but WMO 
Code 75, This is discussed in the two following paragraphs. 

Both references (6) and (16) agree with what many 
oceanographers seem to feel in that at-sea wave height 
observers tend to overestimate wave height. Reference 
(16) even goes so far as to instruct the observer to 
record the average of the estimated heights of the 
larger, well-formed waves. 

Consequently, it seems to make a great deal of 
sense to specify the assumption that observed wave height 
equals significant wave helebt, Hy /3, where Hy/3 i8 the 
average height of the 1/3 highest waves in a seaway. This 
would seem to bridge the gap between the designers who . 
have been prone to use such scales as the one shown in: Figure 
26 because of its mathematically desirable use of Hy /s and 
the ships at sea who have been making observatims based 
on WMO Code 75, 

It might be pointed out that the H, /s Values of 
6, ll, 17, 25 and 37 feet used in this thesis for pre- 
dictions of stern acceleration were chosen as being 
representative of WMO Code 75 Sea States 4, 5, 6, 7 and 


8 respectively. 
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{ WIND AND SEA SCALE FOR FULLY ARISEN SEA 
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